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I. INTRODUCTION 


We introduce ARES, a prototypic system for real-time operational and 
tactical decision support. ARES is designed to quickly and effectively 
help respond to complex emergent problems in disaster relief, in the 
operational art and tactics of warfare, and in related multiperiod, 
large-scale employment of heterogeneous, substitutable resources 
restricted in availability and demand over time, over geography, and by 
organizational limitations. 

Although a great deal of work has been done in strategic modelling in 
many contexts, there is relatively little available modelling help beyond 
simple thumb rules for the time-pressed (operational or tactical) decision 
maker to translate strategic goals into logistically constrained 
operational and tactical plans (and the issues are different). The 
luxuries of hypothetical additional resources and the time to analyze 
their employment are just not available in the operational and tactical 
domains: operational and tactical decisions must be made quickly, and 
usually involve employing only resources actually available to perform 
whatever mission is at hand. 

The history of assignment and allocation models for planning emergency 
logistics extends back to some of the earliest work in linear programming, 
game theory, and their economic interpretation. We cite only a few of the 
references in this large body of literature. The seminal works by Dantzig 
and by Koopmans (both found in Koopmans 1951) are explicitly motivated by 
large-scale logistics problems. Karchere and Hoeber (1953) give early 


direction on the use of newly developed optimization technology in weapon 


system planning and allocation, discussing substitutability of resources 
and choice of suitable objective functions. Geisler (1959) reports RAND’s 
first use of man-machine simulation of logistics support activities. 
Chaiken and Larson (1972) state some basic issues in logistic location and 
task assignment for emergency service vehicles: how many units should 
there be, where should they be located, who should they serve, and how can 
they be relocated to substitute for units not available. Kaplan (1973) 
redeploys divisible resources with linear programming. Fitzsimmons (1973) 
states a nonlinear response time model and uses pattern search to locate 
units well and allocate workload equitably. Swoveland, Uyeno, Vertinsky, 
and Vickson (1973) employ simulation and human interaction to set up a 
unit location problem as a quadratic assignment model which is solved with 
an elegant heuristic. Bracken and McGill (1974) formulate strategic force 
planning models as two-sided games solved with nonlinear programming. 
Bracken, Falk, and Karr (1975) apply multiperiod, two-person zero-sum 
games formulated to develop strategies for unit sortie allocations. 
Finally, Kolesar and Walker (1974) develop a multi-stage solution approach 
to unit and task assignment using set covering and transportation-like 
integer linear programs which are used in real time by applying 
heuristics. 

Named for the Greek god of war, ARES is a proof prototype of a 
real-time decision support system. It employs optimization and simulation 
to capture and exploit a high degree of realism without demanding 
unreasonable amounts of data, or locking the decision maker out of the 
decision process. The intent is to provide quick credible advice with 


good global perspective at a cost no greater than the relatively myopic 


decision methods now widely used. 

ARES accommodates enough detail to support realistic decisions, but 
not so much as to render the process useless. For the intended 
applications, the particular missions to be performed will not likely be 
known much in advance, but the generic types of missions are known and can 
be planned. ARES uses a taxonomy of prepared standardized icons for data 
describing possible missions. The idea is to help the decision maker 
quickly assemble a data scenario closely resembling the proximate 
situation from a menu of these standard icons. 

We characterize the mission at hand as a set of geographically 
dispersed tasks, each composed of partially-ordered sub-tasks requiring 
over time varying amounts of different resources. Organizational units, 
also geographically dispersed and each possessing a different endowment of 
resources, are to be assigned responsibility for the tasks. 
Responsibility for each task rests with only one unit at any given time. 

ARES consists of several models coordinated by a time interval 
simulator which also scales and manipulates scenario data in a fashion 
transparent to the decision maker. Two integer linear programs, a linear 
program, a georeference system, a mobility system, a decision-maker 
simulator, and extensive user interface and user override and control 
facilities complete the program suite. The models in ARES all use a 
standard data interface visible to the decision maker; this invites 
expansion with new models and features. : 

A scenario is created from the attributes of available units and task 
attributes derived in large part from standard cataloged data icons for 


similar tasks. A georeference system is accomodated to generate distance 


costs and delays in relocating and operating units. The decision maker 
may preview the scenario and modify data or manually pre-assign tasks and 
units as he sees fit. 

Operational assignment of tasks to units uses one of two integer 
programming models (IP, ) or (IP). Good task aggregations for the unit 
assigned reduce unit relocation costs and match unit resource endowments 
with aggregated task resource requirements. Logistical considerations are 
paramount at this stage. 

The decision maker can review the operational assignments, modify them 
manually, or reject them outright and restate the conditions for the 
original operational assignment scenario. An acceptible set of 
operational assignments is passed forward to a tactical model. 

Tactical allocation of the resources of each unit to the requirements 
of its assigned tasks uses a linear programming model (GN). Substitutions 
among resources are permitted, although at reduced efficiencies in 
completing the tasks. Allocations recognize task priorities and the 
logistical effects of geographic proximity. In addition, unit efficiency 
in performing a particular task improves over time, and the sequence 
within tasks of resource requirements is considered. The result isa 
complete plan for each unit, showing what resources are to be used to 
fulfill each task requirement, and the efficiency with which operations 
are expected to be carried out. The allocation also determines which 
requirements will not be met in situations which overtax units. 

Finally, the decision maker is presented with a complete solution, 
which he can accept, or modify, or reject outright and repeat. Regardless 


of his action, ARES is designed to lend quick insight. The decision maker 
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uses his own judgement concerning non-quantified factors, and he should 


gain a deeper understanding of the situation at hand from ARES. 


i 


II. OPERATIONAL ASSIGNMENT MODEL (IP) 


This integer program finds good aggregate assignments of tasks to 
units without explicit consideration of unit movement. 


Index Use 


i Tasks 
j Resources 
k Units 


Given Data 
d. Distance cost from unit k to task i 


Minimum, maximum resource j requirements of task i 


=1j ~ 23 

ae a Minimum, maximum resource j employable by unit k 

Zitat Za Penalties for violating minimum, maximum resource limits 
P; Priority of task i (00) 

uss re Penalties for not assigning, double-assigning task i 
Sas Substitution efficiency of resource j (0) 
ae Sequence of resource j requirement within task i 
ne Consumption by task i of resource j from unit k 


Decision Variables 
Xp Binary variable for assigning task i to unit k 


Formulation 


MIN 23 d. Xx 


fe ae tuk | 
Sate ; pe ls!) (eau for all i (1) (GUB) 


Oo nae = . 
As ik Xi. = (eee (ZZ) for all j,k (2) 


10, 1} forall i.k @ (IP) 


ik 
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The notation = (ez) indicates lower and upper ranges (ia) on 
row functional values with corresponding respective linear penalties per 
unit of violation (C425 i.e., this is a goal program with linear 
penalties, an elastic integer program (Brown and Graves 1975). 

Constraints (1) encourage assignment of each task to some unit, and forma 
generalized upper bound (GUB) row set (Dantzig and Van Slyke 1967). 
Constraints (2) express the goodness of fit of task assignments with 
employable unit resources. Constraints (3) preclude fractional assignment 
of tasks to units. 


Consumption by task i of resource j from unit k is defined: 


<2 
= —In f.. + (p,-1)/10 + d., fo, + t, 
Me = ia i ik ~k ik , (1) 


where O71. is the speed of advance of a unit and ty is the number of 
periods that unit k has already been assigned task i. The rationale for 
the particular consumption function (1) amplifies the resource requirement 


My to account for the state of resource readiness fay the task priority 


P; (making less important tasks appear more expensive), the logistic 


proximity of unit k and task i, d 1/1: and learning curve effect as a 


function of time since assignment, t The data are scaled so that (1) 


ik” 
is in conformity with policy guidance or the judgement of the decision 
maker. Alternate consumption functions may appeal in other situations. 

The distance costs d. 


and penalties ue, Uy and Zz. are expressed 


k jk’ * jk 
in commensurate units and deserve some thought by the modeller. For 


instance, Za may be interpreted as how much additional distance cost 


should be incurred before considering overtaxing maximum resource 
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employment a for unit k; this is a direct expression of logistical 
efficiency. For simplicity in our tests, distance costs di. are scaled by 


a policy parameter, ar and Za are part of the input script, u. is 


defined as 100/p., and u, equals 100. 
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III. OPERATIONAL ASSIGNMENT MODEL (IP, ) 
The purpose of this integer program is to find good movements of units 
to locations from which they will be assigned good aggregate groups of 


tasks to perform. 


Index Use 
i Tasks 
j Resources 
k Units 
l Locations (assumed here to be collocated with tasks) 
Given Data 
diz Distance cost from unit k to location 1 
B51 Distance cost from task i to location 1 
a Gross resource requirement j of task i performed 
from location 1 
ii Net resource availability j of unit k located at 1 


Decision Variables 


Z Binary variable for assigning task i to location l 


ry 


x Binary variable for moving unit k to location l 


lk 
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Formulation 


MIN >> 8.) 73117 2d x 


«1 iy | kl [ke ell 
S Z. = (1,1):(u,.u,) for all eG@micies 
1 wl el 
» He = (1,1);(m,m) for all k (2) (GURR 
] 
> Koo) =e (Oe) aC mern) for all 7 i 
7 lk 
o e 
=2y nt : Xi, = (O,1);(m,m) for all i,1 (4) 
= Tal 741 + aie ee (G70), ) for a2ael,j 
Z +4 = One for all 1.15 
Xi, = {0,1} for all lke (IP, ) 


(IP, ) uses the notation of (IP). Constraints (1) encourage assignment 
of each task to some location. Constraints (2) allow movement of each 
unit to some location. (A GUB row set is formed by constraints (1) and 
(2).) Constraints (3) attempt to restrict assignments so at most one unit 
is moved to any particular location, (4) require that a unit be moved to 
any location to which a task is assigned, and (5) attempt to match for 
each location and each resource an aggregate assignment of tasks which 
have gross resource requirements about equal to the net resource 
availability of the unit moved to that location to perform the tasks. 
Constraints (6) and (7) preclude fractional location of tasks and units. 


Gross resource requirement Tail represents the resource j estimated to 


be required at location 1 in order that task i actually receive Tia 


ae = —In f .20( pe) 10 ee 
Aol a Lae e i3 i i) se (2) 


where o expresses the logistic radius of influence from any location; we 
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have used o = 100. The gross resource requirement (2) amplifies the 
resource requirement tae in the same fashion as (1). 


Net resource availability a. represents the amount of resource j 


jik 


which unit k can deliver from its endowment 2. forward to location l. 


Unit k may be moving toward location 1 while supplying this net resource. 


ay ay =d 70. 
ar ike = Gabe fi, (a1. + (1 a1.) e 1k ‘k), (3) 


where G1, is the fraction of time which the unit will spend at its 


destination location, and o 


k is the speed of advance. 


The distance costs diy. and Bi. and the penalties us toa m, b and b 
all render the same objective function units. In our work, m = 100, and 
u, and a are defined as in (IP). The penalties for assigning too little 


(or too much) resource j to location 1 are b (or b). We have used b=0.1 


and b = 0.01. 
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IV. TACTICAL ALLOCATION MODEL (GN 


) 


This linear program allocates resources to the tasks assigned to unit 


Ke 


Index Use 


J 
W 


Given Data 


e. 
iwj 


Tasks 


Resources 


Work (resources required by tasks assigned to unit k) 


Minimum, maximum work requirements w of assigned task i 


Penalties for violating minimum, maximum work requirements 


Minimum, maximum resource j employable by unit k 


Penalties for violating minimum, maximum resource limits 


Priority of task i 


Substitution efficiency of resource j for work requirement w 


(0) 


Sequence of work requirement w in task i (D0) 


Efficiency of resource j used for work w on task i 


Decision Variables 


Viwj 


Formulation 


MAX 22> 
st iwj 
> 

J 

23 

iw 


Allocation of resource 


“iwj Yiwj 
re) 
Yiwj 
“iwj Yiwj > 
Yiwj 2 


j to task i resulting in work w 


(iy? iw) 


, ee 
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for all i,w (1) (GUB) 
for alia (2) 


for all i,w,j (3) (GN, ) 


(GN) uses the notation of (IP). However, the dimensions of (GN, ) 
discriminate between resources consumed, j, and work completed, w, 
explicitly representing substitutability of resources. Constraints (1) 
encourage allocation of sufficient work resources, while constraints (2) 
indicate the desired mix of employable unit resources. Constraints (3) 
require non-negative resource allocations. (GN) is an elastic 
generalized network (Brown and McBride 1984). 


Efficiency of resource j used for work w on task i is defined 


+ = 
23 peur ta + (p,-1)/10 + S 54/10 + dif. (4) 


Saar 
iwj 


+ : , 
where Si, = max {O, s. - t}, t is the last time period of this 


1W 


allocation, and 1, 


employs the readiness and substitutability of resources via # oct S 


is the unit speed of advance. The efficiency (4) 
+ 
1W 


reduces efficiency if the work w should not be started until period Sy 


WS, 


If model (IP) has been used for operational assignment, 


ik = max {O, di. 7 01/2}. (5) 


If unit k is to be advanced toward, or to, location 1 by model (IP, ). 


di. = Macc Or dix - 01/2} + .)- (6) 
These distance costs dix in (5) or (6), and penalties Diw? Tiw? a and 
Zo are all intended to yield the same objective function units. For our 
_ + - _ 
tests, Ca 100/p.s._. and diy = 100. 
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V. QONSIDERATION OF LOGISTICS 


The efficiency with which a unit completes a task depends heavily upon 
logistical considerations. If a unit is remote from a task, or must be 
moved, its efficiency suffers. Figure 1 shows an idealized situation with 


unit k, task i, and location 1. 


(ease se 
t I 
| | 
irae 

SS el 

| | 

\ | 

L.-J omer 9 
(Saeaiont 
1 1 
! 1 
| 





Figure 1: Idealized Geographic Logistical Scenario 


Model (IP) assigns tasks to units relying exclusively upon d. (IP 


kK" i) 
moves units to new locations and assigns tasks to be performed from these 
new unit locations. (IP, ) recognizes d,, and g.,. The distances d., and 
L lk tl ik 

dip are surrogates for logistical costs of assignment during the ensuing 
time period. Clearly, (IP, ) is more appropriate for situations in which 
unit movements are expected, (IP) when they are not. (IP, ) provides the 
decision maker with a better opening gambit than does (IP) if the scenario 
involves significant initial redeployment of units. 


Tactical allocation models (GN) are given unit and task assignments 


and planned unit movements. Therefore, (GN) can allocate resources using 


PA) 


any logistic efficiency function of assigned distances, and of other 
attributes induced only from assignment such as weather effects, speed of 
unit movement, etc. (GN) can also substitute resources at somewhat 
reduced efficiency as well as prioritizing their immediate application. 
Given a fairly reasonable operational assignment, (GN) provides a 
high-resolution work plan with rich logistic detail and good face 


validity. 


ae 
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AN EXAMPLE SCENARIO 


We demonstrate ARES with an example constructed for Engineering 


Battalions of the Hellenic Army. 


repairing major damage to public works following an earthquake. 
purposes, there are 14 units, 


Figure 2 shows the units and tasks from the ARES input script. 


The mission scenario involves 20 tasks 


For our 


each endowed with some of 25 resources. 


In the 


United States, the Department of the Army defines unit types in (1976), 


and task standards in (1973a). 
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: Units and Tasks of Example 


Wy ILABEL 
1 ADMIN BUILDING 441051 
2 ADMIN BUILDING 44105] 
3 ADMIN BUILDING 44110] 
4 HOSPITAL 100 BED GHOI1]]) 
5 HOSPITAL 200 BED GHOCTI 
6 HOSPITAL 100 BED GHO1ZI 
7 HOSPITAL 100 BED GHOI31 
B RAILROAD BRIDGE B61643 
9 RAILROAD SRIDGE 861512 
10 ROAD BRIDGE $0° 8sSc¢I0I 
11 ROAD BRIDGE 100° 884109 
es ROAD BRICDSE 70° BS4106 
13 ROAD 2.5 MILES 853120 
14 ROAD 6.7 MILES BSeslew 
15 ROAD §.5 MILES BSi128 
16 ROAD 6.8 MILES 853126 
17 ROAD 6.0 MILES BSZI20 
18 WATER TANK-DIST-SUP NOI 
19 WATER TANK-DIST-SUP NO2 
a0 WATER TANK-DIST-SUP NO3 
Figure 2 


The georeference coordinates of units and tasks are given in Figure l 


for the situation depicted in Figures 3, 4 and 5. 
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Figure 3: Initial Geographic Locations of Units. (Coordinates 
displayed are a georeference in common with the following 
figures. ) 
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Figure 4: Earthquake Epicenter. 
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Figure 5: Geographic Locations of Tasks. Geographic locations of 
damaged public works and earthquake epicenter are shown. 
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A georeference system is used to generate coordinate-to-coordinate 
distance costs, which appear in the ARES input script. 

The resource requirements for Task 1 ("ADMIN. BUILDING AA1051"), a 
disaster relief facility, are given in Figure 6 and the input script. 
Resource requirements such as these are available in standard engineering 
reference manuals for a wide variety of task types (for instance, see 
unclassified sources from the United States Department of the Army (1973a, 
1973b, and 1973c). We envision a taxonomy of standardized task data icons 
from which a particular set of requirements can be very quickly extracted 
and assembled for a scenario. The size of our resource requirements data 
base is modest but the resulting accuracy and level of detail are quite 
good. Better yet, data mobilization from a menu of such icons can be 


completed in minutes. 


RESOURCE LABELS AND TASK 1 REQUIREMENTS 


nN WwW 


ROTARY TILLER OPER, 
SCRAPER OPER. 


RR iLABEL { HOURS 
1 ENGIN=-P LON-APREN-HLPER 6648 
2 SURVEYOR 70 
3 CARPENTER 7557 
% ELECTRICIAN 960 
S PLUMBER 1760 
6 MASON 1600 
? STRUCTURE SPECIAL. 0 
8 HEAT-VENTILAT SPECIAL. 200 
9 WELDER 0 

10 PIPELINE 

11 CRANE-SHOVEL OPER. 200 

I LOACER OPER. z7S0 

ie DOZER PER. 300 

16 COMPRESSOR OPER. 0 

1S DUMP TRUCK OPER. 600 

16 CONCRETE MACHINE OPER. 0 

1? GRADER OPER. 5298 

18 CRUSHER OPER. 0 

19 DITCH MACHINE OPER. ~ 500 

20 ASPHALT SPECIAL. 

all POWER ROLLER OPER. 

22 WATER DISTRIBUT. OPER. 

2 


0 
G 
0 
POWER BOAT OPER. 0 
0 
0 


rd 
wt 


Figure 6: Resource Requirements of Task 1 
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The resources employable by Unit 1 ("1ST COMBAT BN"), a combat 
engineering battalion, are given in Figure 7 and in the input script. 
These resource endowments are in line with those given by the United 
States Department of the Army (1973a) with conversion to man hours from 


flor). 


RESOURCE LABELS AND UNIT 1 AVAILABILITIES 


RR ILA8EL i 'MIN iMAX IMiN PEN IMAX PEN 
e SURVEYOR 405 <S0 10 10 
3 CARPENTER Tabs 1350 10 10 
4 ELECTRICIAN =035 ceo 10 10 
5 PLUMSER 1418 LS75 10 10 
6 MASON Pees 1250 19 10 
a Fai Se ol): | a ee os ee 0 0 10 10 
8 HEAT-VENTILAT SPEC.AL. 0 0 10 10 
S WELDER 187 208 10 10 

10 PIPELIUNME 0 0 10 10 

jl CRANE-SHOVEL OPER. rs 1250 10 10 

le LOADER OPER. oles 6850 10 10 

LS iCeER OPER: 4050 6500 10 10 

14 COMPREZTICR OPER. 1013 Li2s 10 10 
15 DUM TRUCK OPER. 10935 12150 10 10 

16 CONCRETE MACHINE OPER. 203 ies 10 10 

Ws GRACER OPER. 1620 1800 10 10 

18 CRUSHER OPER. 0 0 10 10 

Ms) DITCH MACHINE OPER. 0 0 10 10 

20 ASPHALT SPECIAL. 0 0 10 10 

21 POWER ROLLER OPER. 0 0 10 10 

a WATER DIZTRIBUT. OPER. 0 0 10 10 

a5 PCKER BOAT OPER. 0 0 10 10 

26 ROTARY TILLER OPER. 0 0 10 10 

<5 SCRAPER OPER. 0 0 10 10 


Figure 7: Resource Endowment of Unit 1. 


The input script also includes for each task the sequence of resource 
requirements expressed as the first period when the resource is best 
applied, and for each resource its substitution efficiency for other 


resources. 


The scenario data constitutes about 1,000 records. However, these 
records derive from the unit, task, and resource definitions which are 


modest in number. 
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VII. DESIGN AND IMPLEMENTATION 


ARES is intended to help the decision maker, not to replace him. 
Figure 8 shows the functional structure of ARES. The design is biased 
toward interactive use with review and intervention options at each stage 


of operational assignment and tactical allocation. 
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cNITTALIZE: 
NEXT_PERIOD: 


OP_ASSIGN: 


REVIEW_IP: 


TAC_ALLOC: 


UNIT-K: 


NEW_SCRIPT: 


REVIEW_PERIOD: 


Define NEW_SCRIPT 

Redefine NEW_SCRIPT as OLD_SCRIPT 

Select Model (IP, ) or (IP) 

Read OLD_SCRIPT 

Generate and Solve (IP, ) or (IP) 

Record task and unit assignments on ASSIGN_FILE 
Option to review assignments in ASSIGN_FILE 
either stop, 

or edit OLD_SCRIPT and GOTO OP_ASSIGN, 

or edit OLD_SCRIPT and/or ASSIGN_FILE and continue 
Read OLD_SCRIPT and store as SCRIPT 

Read ASSIGN_FILE and update SCRIPT assignments 
Select (GN, .), Generate and Solve 

Update SCRIPT resource requirements for work completed 
For next unit k REPEAT UNIT-K 

Update SCRIPT unit locations and distance costs 
Write SCRIPT as NEW_SCRIPT 

Option to review results 

either stop, 

or edit OLD_SCRIPT and/or ASSIGN_FILE 

and GOTO OP_ASSIGN 


or edit NEW_SCRIPT and GOTO NEXT_PERIOD 


Figure 8: ARES Functional Specification 
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ARES is implemented in FORTRAN H(Extended) and executes on an IBM 3033 
AP computer using the VM/CMS operating system. Input scripts are read 
from files which may be viewed and modified with a full-screen editor such 
as XEDIT. (Software copyrights IBM Corporation. ) 


ARES uses the X-SYSTEM (Brown and Graves 1975) to solve (IP (IP) 


1) 
and (GN, ) in real time. For each problem instance, problem generators 
directly convert input script data into an internal representation, the 
solver is invoked, and the solution is provided to a report writing 
program. ARES consists of a set of open subroutines and is executed with 
whatever preview, review, or other external interference is deemed 
desirable. 

We envision cyclic use and review at varying levels of detail asa 
mission progresses over time. Accordingly, input scripts include the 
beginning period and number of periods in the ensuing time interval, which 
intrinsically scales time-dependent input data to the desired level of 
aggregation. We have tested ARES manually and by replacing the decision 
maker with a simulation which performs “judgement review" of successive 
solutions over time. This permits totally automatic evaluation of 
complete mission scenarios, and avoids tedius manual effort in our 
research. (A single time interval may generate 15, or 20 thousand lines 
of solution detail at the scale of our example scenario. ) 

The update of unit coordinate locations and distance costs is a simple 
surrogate for a more realistic and complicated georeference and mobility 
system. ARES estimates the direction and speed of advance of each unit 
during the time interval and relocates the unit. Then the distance costs 


are adjusted. If operating areas are known sufficiently in advance to 
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permit preparation of detailed georeference and mobility systems, ARES can 
accomodate the increased level of detail in real time (e.g., Brown, Ellis, 
Graves, and Ronen 1987). The update can also be used to degrade, or to 
amplify unit resource endowments and effectiveness to modify task resource 
requirements, or to change any other data artifact, providing a rich 


modelling arena. 


on 


VIII. SCENARIOSRESTETS 


ARES has been used in simulation mode to completely plan mission 


scenarios from start to finish. For the earthquake scenario, Figure 9 
shows the initial operational assignments of (IP, ). 


Tia 
ly 


ce 


Ce 








Figure 9: [Initial Operational Assignments of Units. Directional 
vectors show the straight-line path and relative speed of 


advance 1: 
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Figure 10 depicts the arrival of units to their initially assigned 


locations. 





Figure 10: Initial Operational Assignments of Units to Locations. 
Arrows show straight-line path of advance toward assigned 
locations. 


Without intervention by the decision maker, the scenario completed 
itself in 7 weekly intervals, requiring less than 2 minutes ina 1.2 
megabyte memory region. 

Face validity of the simulation solution is excellent. No manual 
intervention has been found to improve the solution. In fact, many manual 


attempts to coerce better assignments resulted in startling degradations. 
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The application of available resources, with allowable substitutions, 
is shown in Figure 11 for the 7 single-period time intervals required to 


complete the earthquake scenario. 
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Figure 11: Resource Requirements and Work Completed. Each row 
represents a resource requirement over time-interval columns. 
The white bars depict resource requirements by time 
interval; the black bars show the relative fulfillment of the 
requirements. Broken bars are out of scale. From each time 
interval to the next the requirements are reduced by the work 
completed and amplified by new sequence-dependent 
requirements. In this scenario, ¢ weekly time intervals are 
required to complete all tasks. 


34 


IX. COMPUTATIONAL EXPERIENCE 


Extensive computational experience reveals that the operational 
assignment models (IP) and especially (IP, ) are most difficult to solve at 
the beginning of a scenario, and get progressively easy in later time 
intervals. The size of these models varies with the number of mandated 
assignments, impossible assignments, and the non-zero density of resource 
availabilities and remaining requirements. (IP) typically has about 340 
constraints, 268 binary variables, and 6,200 non-zero consumption 
coefficients. The linear program continuous relaxation can be generated 
and solved in about 5 seconds, and an optimal binary solution is achieved 
in another second, or so. 

(IP, ) has about 1,000 constraints, 645 binary variables, and 8,000 
rather unwieldy non-zero gross resource requirement and net resource 
availability coefficients. 


The linear program continuous relaxation of (IP, ) proved impossible to 


1) 
solve by direct assault. Prior work by Brown and Graves for Bausch (1982) 
on large-scale set partitioning problems and later refinements by Brown, 
Graves and Ronen (1987) suggested an alternate means of attack: a problem 
cascade. 

Briefly, the rows of constraints and columns of variables are 
lexicographically sorted to place short rows first accompanied by other 
rows and columns with intersecting non-zero coefficients, and longer rows 
later with their own intersecting rows and columns. 


The problem cascade proceeds by activating a set of constraints, 


relaxing all other constraints, and activating a set of variables, fixing 


all other variables to their last-known values. This problem is solved, 
the new values of the active variables recorded, and another problem 
specified in the building problem cascade. The last problem in the 
cascade activates all constraints and variables (precisely the problem 
found intractable above) and solves it by starting with an advanced 
solution recorded from the last-known values of variables solving previous 
problems in the cascade. 

(IP, ) resisted even the problem cascade until a new heuristic cascade 
strategy was adopted which activates the shortest 1/2 of constraints and 
their associated variables, then the shortest 3/4, then 7/8, and so forth 
until the last constraint is added and the problem is solved. Remarkably, 
this approach has been absolutely reliable and robust, while most others 
fail or prove unruly. 

Generation and complete problem cascade solution of the continuous 
relaxation of (IP, ) now requires about 10 seconds. 

An acceptable binary solution to (IP, ) is achieved in another second, 
or two. 

We do not routinely seek optimal binary solutions to (IP, ). which we 
refer to as "perfect misfits”. The gross resource requirements and net 
resource availabilities in (IP, ) are rough logistic estimates, calibrated 
by actual field experience but ultimately just approximate target 
performance levels. For interesting operational assignments (i.e., early 
in the scenario) there are simply no feasible solutions; the goal is to 
guess where to send units so that they can peremptorily cope with the 
mission at hand with maximal effectiveness. Accordingly, we accept in 


practice binary solutions which may be as much as 25% greater than an 
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optimal lower bound in total value, including constraint violation 
penalties. Experimentally, we have determined at additional computational 
cost (as much as 5 minutes per trial) that these binary solutions are 
actually almost always within a few percent of the true optimum. 

A decision maker can help ARES with its operational assignments or 
completely specify a solution with manual assignment features. Our 
experience suggests that the decision maker can express some 
non-quantifiable guidance in this fashion, but can not hope to apply a 
remotely competitive global perspective. Manual competition with ARES 
reveals that model computation effort is amply justified by the quality of 
operational assignments achieved. The operational assignment models, 
especially (IP, ), produce solutions no decision maker is likely to 
discover. Some of these solutions have yielded remarkable insights. The 
initial operational commitment of units is arduous and crucial to mission 
success. | (IP, ) is worth the computational investment. 

By contrast, the tactical allocation models (GN) are easy to solve 
even in the cases where heroic substitution of resources are required. 
The size of each (GN) varies with the number of tasks assigned to the 
unit, and the non-zero densities of resource availabilities, remaining 
requirements, and allowable substitutions. For our scenario, a typical 
instance of (GN) has about 70 constraints and 1,190 variables, and is 
generated and solved in less than 0.1 second. Stress tests with 525 


constraints and 12,500 variables require less than a second. 
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X. DISCUSSION AND CONCLUSION 


The subtlety of operational assignment has surprised us, as has the 
ease of detailed tactical allocation. Operational assignments are 
delicate decisions, and the success of entire missions appears to be very 
sensitive to minute details---precisely the considerations a hard-pressed 
decision maker would likely overlook in haste. 

Extensive mechanisms have been provided in ARES to encourage manual 
review and coercion of solutions. However, there have been very few cases 
in which such guidance improved solutions and many instances in which 
minor manual adjustments of operational assignments inflicted great 
disruption. For example, some operational assignments of (IP, ) 
“cross-locate” units in the sense that a pair of units will each be 
collocated with a task assigned to the other. This superficial blemish 
can easily be masked by manual intervention or by automated solution 
editing. Surprisingly, the removal of cross-locations frequently 
increases the logistic cost of the solution: there is a very delicate 
balance of logistic support of task cohorts assigned to specialized units. 
Cross-location can actually make a great deal of sense in practice. 

Manual intervention can work well in cases inviting human judgement. 
For instance, nearly completed tasks or tasks which have been in progress 
for long intervals can enjoy efficiencies not apparent to our models. The 
decision maker can easily declare tasks completed when minor requirements 
remain, or when it is clear that the models are unduly influenced by a 


minor requirement. 


Operational assignments can be restricted so that units are not moved 
from their initial new locations until the work in their logistic 
influence has been completed. Surprisingly, this restriction is rarely 
needed in practice, and in those cases in which multiple relocations are 
indicated great efficiencies accrue to the mission as a whole. We view 
this insight as a strong validation of the modelling philosophy underlying 
ARES. 

Fortuitous design decisions to separate operational assignment and 
tactical allocation models, to decompose time intervals, and to couple the 
resulting restricted components with simulation and human intervention 
options have yielded more than the intended benefits. Our original 
motives were to capture as much reality as possible while still rendering 
models capable of quick, responsive solution. 

The decomposed design also naturally accomodates features which are 
otherwise difficult to provide. For instance, partial orderings within 
tasks can be introduced. Also, discussions with Professor Wayne Hughes 
have suggested the technical feasibility of campaign analysis, two-sided 
gaming, and force-on-force applications of ARES. In these contexts, the 
coupling with simulation enhances our capabilities enormously. 

ARES was originally designed for use on an IBM-PC. There is no 
compelling reason not to use this microcomputer, but we encountered a few 
practical limitations. An arbitrary 640 kilobyte memory region limitation 
and crippling errors in the FORTRAN compilers available for the IBM-PC 
present artificial conversion costs which we are not willing to bear. 

When these unfortunate shortcomings of the IBM-PC are repaired, conversion 


might be reconsidered. Calibration tests project solution times on the 
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order of 2 minutes per time interval on IBM-PC/AT with a math co-p 


_ 


and internal clock speed of 8 megahertz. - 
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